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I.  INTRODUCTION 


For  many  years  computer  specialists  have  been  trying  to  model 
physical  phenomena  by  using  hydrodynamic  computer  codes.  One  such 
phenomenon  is  the  shaped  charge  liner  collapse.  This  report  is  a 
comparison  between  ~a  hydrodynamic  computer  code  and  experimental 
results  of  this  phenomenon. 

The  hydrodynamic  computer  code  used  in  this  comparison  is  the 
BRLSC  (Ballistic  Research  Laboratory  Shaped  Charge)  code1.  This  code, 

a modified  version  of  the  HELP  code2,  is  Eulerian  formulated;  that  is, 
a computational  grid  or  coordinate  system  is  fixed  in  the  region  of 
interest,  and  the  flow  variables  are  calculated  for  each  cell  in  the 
grid  as  the  specified  material  moves  through  the  cells.  This  code  was 
developed  specifically  for  the  solution  of  the  shaped  charge  problem 
which  is  two-dimensional,  unsteady,  compressible  flow.  The  BRLSC  code 
was  developed  by  System,  Science,  and  Software,  Inc.  under  contract 
with  the  Ballistic  Research  Laboratory  (BRL) . 

The  experimental  results  are  from  two  shaped  charge  systems.  The 
first  system,  refered  to  as  the  "43mm  shaped  charge  or  43mm  cone  di- 
ameter", is  a hollow,  conical,  copper  liner,  driven  by  a right  circular 
cylinder  of  bare  composition  B explosive.  The  second  system  refered  to 
as  the  "105mm  shaped  charge"  is  also  a hollow,  conical  copper  liner, 
driven  by  an  unconfined,  right  circular  cylinder  of  composition  B 
explosive. 

The  comparisons  will  consist  of  various  experimental  results 
(such  as  early-time  flash  radiographs  of  the  liner  collapse,  jet  tip 
velocity,  jet  and  collapse  velocity  distributions,  and  jet  density 
and  temperature  measurements) , and  their  relationship  to  theoretical 
computed  data  presented  in  graphical  form.  Another  series  of 
comparisons  will  be  presented  from  a calculated  parametric  study  con- 
cerned with  the  variations  in  the  equations  of  state  of  both  liner  and 
explosive.  A detailed  analysis  of  the  calculated  results  will  also 
be  presented. 

The  objective  of  this  report  is  to  present  results  which  will  show 
the  warhead  designer  the  accuracy  with  which  the  BRLSC  code  can  predict 
the  liner  collapse  and  jet  formation  process. 


1A/.  L.  Sittings,  "BRLSC:  An  Advanced  Eulerian  Code  for  Predicting 
Shaped  Charges,  Volume  I, " Ballistic  Research  Laboratories  Contract 
Report  No.  279 , December  1975.  (AD  M023962). 

2L.  J.  Hageman  and  J.  M.  Walsh , "HELP,  A Multi-material  Eulerian 
Program  for  Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two 
Space  Dimensions  and  Time,  Volume  I,  "Ballistic  Research  Laboratories 
Contract  Report  No.  39,  May  1971.  (AD  #726459). 
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II.  BRLSC  CODE  DESCRIPTION 


BRLSC  is  a FORTRAN  code  which  has  evolved  from  the  OIL3  hydrodynamic 
code  and  developed  over  a 15  year  period. 

The  system  of  partial  differential  equations  which  is  solved  by  the 
code  consists  of  the  balance  equations  of  mass,  momentum,  and  energy 
written  with  respect  to  a coordinate  system  fixed  in  space  (Eulerian) . 
These  equations  along  with  an  equation  of  state  and  a stress  - strain 
relationship  are  sufficient  for  the  numerical  solution  of  the  detonation 
of  high  explosive,  shaped  charge  liner  collapse,  and  jet  formation. 

BRLSC  code  is  time' dependent  in  two  space  dimensions  with  an  option 
for  either  plane  (x,y)  coordinates  or  axisymmetric  cylindrical  (r,z) 
coordinates.  The  time  increment  is  based  on  the  Courant  stability 
condition  and  the  maximum  sound  speed,  plus  the  particle  velocity  in 
the  grid.  The  resulting  numerical  technique  is  an  explicit,  conditionally 
stable,  finite-difference  system  which  is  used  to  integrate  the  governing 
equations. 

Solutions  to  these  shaped  charge  problems  involve  the  motion  of  a 
metal  liner,  which  is  accelerated  by  the  pressure  loading  from  a 
detonating  explosive.  The  gas-metal  interface  is  separated 
by  lines  joining  massless  tracer  particles.  These  particles  move 
across  the  grid  and  locate  the  interfaces  and  free  surfaces  of  each 
material.  Consequently,  a Lagrangian  effect  is  introduced  into  the 
treatment  of  moving  surfaces. 

The  BRLSC  code,  using  the  Von  Mises  yield  condition,  has  as  an 
option  the  treatment  of  material  as  being  elastic  - perfectly  plastic. 

But  for  the  results  reported  herein,  the  liner  material  is  treated  as 
being  purely  hydrodynamic. 

Plotting  programs  have  been  developed  by  the  author  at  BRL 
to  display  the  voluminous  data  generated  during  a computer  run. 

In  those  figures,  as  indicated,  the  tracer  particles  were  used  to 
outline  the  interfaces  and  free  surfaces  of  the  different  materials. 

Since  the  problems  reported  here  are  axisymmetric,  only  half  of  the 
various  fields  are  presented. 


3V.  E.  Johnson , "OIL3  A Continous  Tivo-DimensionaZ  EuZerian  Hydro- 
dynamic Code"3  GuZf  GeneraZ  Atomic , G AMD-5580  Revised , January 
1965. 
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III.  COMPARISON  BETWEEN  BRLSC  CODE  RESULTS  AND  43MM  SHAPED  CHARGE 
EARLY-TIME  FLASH  RADIOGRAPHS 

Figure  1 is  a line  drawing  of  the  43mm  shaped  charge,  that  was 
used  in  this  comparison.  Figure  2 shows  the  early-time-flash  radio- 
graphs of  the  liner  collapse  of  this  charge.  These  radiographs  were 
made  by  Boyd  Taylor  of  the  Ballistics  Research  Laboratory  (BRL).4  It 
begins  with  the  arrival  of  the  detonation  wave  at  the  apex  of  the 
cone,  and  continues  at  2.4,  3.8,  5.6,  6.5,  7.6,  14.6  ys.  A final 
exposure  was  made  at  20  ys . The  simulated  shaped-charge  liner 
collapse  sequence  is  depicted  in  Figure  3 at  five  discrete  times. 

These  five  times  agree  with  those  of  the  radiograph  at  3.8,  5.6, 

6.5,  7.6  and  9.6  ys.  However,  the  code  calculations  do  not  include 
all  of  the  time  sequences  from  the  flash  radiographs  shown  in  Figure  2. 
The  first  radiograph,  corresponding  to  2.4  ys,  is  not  presented 
because  the  jet  is  not  clearly  defined  in  the  picture.  Also,  those 
radiographs  whose  times  are  greater  than  9.7  ys  were  excluded  because 
the  computer  run  time  became  excessive.  The  simulated  collapse 
sequence  agrees  very  well  with  those  of  the  flash  radiographs  with 
which  we  can  make  comparisons. 

Figure  4 is  a calculational  plot  of  the  BRLSC  code's  results  of 
the  velocity  difference  between  the  tip  of  the  jet  and  the  tail  of 
the  slug,  as  a function  of  time  after  the  detonation  wave  hits  the 
apex  of  the  cone.  This  difference  is  depicted  by  the  broken  line, 
and  hereafter  will  be  referred  to  as  "relative  jet  tip  velocity." 

This  figure  demonstrates  that  after  approximately  11.5  ys,  the  jet 
tip  reaches  a constant  relative  velocity  of  6.54mm/ys. 

Figure  5 shows  a comparison  between  the  calculations  of  the  BRLSC 
code  in  Figure  4,  and  the  measured  distances  between  the  jet  tip  and 
the  slug  apex  as  a function  of  time,  as  depicted  in  Figure  2. 

The  slope  of  the  line  passing  through  the  flash-radiograph  data 
points  is  a constant,  6.54.  This  slope  is  the  relative  velocity  of 
the  jet  tip  with  respect  to  the  slug's  tail,  6.54mm/ys.  This  is  in 
agreement  with  the  calculations  of  the  BRLSC  code  depicted  in 
Figure  4. 

However,  this  graph  shows  a difference  of  approximately  1 ys  in 
the  movement  of  the  jet  tip,  between  the  BRLSC  code's  results  and  the 
measurements  from  the  radiograph.  This  difference  is  a result  of 
either  a timing  problem  in  the  recording  instruments,  or  numerical 
problems  regarding  equations-of-state  for  the  explosive,  the  metal, 
or  both.  Later  in  this  report,  a study  of  equations-of-state  will 
be  presented. 

^Private  Communiaations  from  B.  Taytor3  Ballistic  Research  Laboratory 3 
Aberdeen  Proving  Grounds , Maryland. 
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DIMENSIONS 
A - 43 
B - 46 
C - 73 
D ■ 1.25  mm 

, • 


Figure  1.  The  line  drawing  of  the  43mm  shaped  charge.  Pertinent 
dimensions  are  shown  in  the  table. 
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SHAPED  CHARGE  LINER  COLLAPSE 

43MM  DIA.,  1.25  MM  COPPER  WALL  CONES 
600 KV,  0.2  USEC.  EXPOSURE  AT  TIMES  INDICATED 


L 


Figure  2.  Experimental  flash  radiographs  of  the  43mm  shaped  charge  liner  collapse 

sequence.  This  figure  is  from  reference  4.  The  unit  of  time  is  a microsecond 
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Figure  3.  BRLSC  code  calculations  depicting  the  43mm,  1.25mm,  copper  wall,  shaped  charge  liner 
collapse  sequence.  Times  are  3.8,  5.6,  6.5,  7.6,  and  9.7us  left  to  right  and  top 
to  bottom  respectively,  with  respect  to  Figure  2. 


Figure  4.  BRLSC  code  results  from  43mm  shaped  charge  calculations. 

The  relative  velocity  of  the  jet  tip  with  respect  to 
tail  (slug)  is  shown  as  the  dashed  line. 
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O EXPERIMENTAL  DATA 
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TIME  (n. s) 

Figure  5.  Comparison  between  the  BRLSC  code  results  and  measured  differences  of  the  jet  tip 
with  respect  to  tail  (slug),  denoted  by  solid  line  taken  from  flash  radiographs 

(Ref.  4) 


Setting  aside  the  small  difference  previously  mentioned,  one  can 
conclude  that  the  simulation  is  accurate  on  the  basis  of  the 
aforementioned  comparisons, 

IV.  DETAILED  ANALYSIS  OF  THE  CALCULATIONS  FROM  THE  43MM  SHAPED 

CHARGE 

The  first  graph  in  a detailed  analysis  of  the  43mm  shaped  charge 
calculations.  Figure  6,  illustrates  the  axial  component  of  velocity 
for  liner  material  along  the  axis  of  symmetry  as  a function  of  axial 
distance.  The  time  intervals  of  8,  10,  and  13ys  are  presented.  The 
so  called  "inverse  velocity  gradient"  in  the  jet  is  illustrated  at  the 
time  intervals  of  8 and  10  ys.  This  phenomena  has  been  observed5 
and  calculated6  earlier  by  using  a one-dimensional  model  of  the 
shaped  charge  problem.  This  inverse  velocity  gradient  in  the  jet  is 
finally  damped  out  by  13ys,  and  the  jet  tip  is  equilibrated  to  a 
velocity  of  approximately  6.8mm/ys. 

The  next  graph  in  the  analysis.  Figure  7,  contains  three  relations 
as  functions  of  the  axial  position,  for  liner  material,  along  the 
axis  of  symmetry.  These  three  relations  are  ratios  of  density,  p/pQi 

specific  internal  energy,  1/ Iq ; and  pressure,  P/P^.  Where  the 

following  are  the  normalizing  units  used  in  the  above  relationships 
as  well  as  those  in  the  remainder  of  this  report  unless  otherwise 
stated: 


a.  pQ  = 8.9  Mg/nf5 

4 

b.  Iq  = 10  joules/g 

c.  P = 1.0  Mbar. 

o 

In  the  vicinity  of  the  jet  tip,  the  density.  Figure  7a,  is  50%  of  the 
liner's  initial  density.  The  sublimation  energy  level,  (1^  on  the 

graph)  Figure  7b,  indicates  that  the  calculated  specific  internal 
energy  of  the  same  region,  i.e.,  jet  tip,  is  six  times  greater.  Also, 
eighty-three  per  cent  of  the  liner  material  along  the  axis  of  symmetry 
has  specific  internal  energy  greater  than  1^.  This  indicates  the  liner 

material  has  undergone  vaporization,  i.e,  the  expanded  hot  state,  when 
p/p  <1  and  I > 1 1 and  transition  state,  when  p/p  < 1 and  I < I<  I*. 

OS  o O s 


SR.  DiPersio  and  J.  Simon 3 "An  experimental  Method  of  Obtaining 
Collapse  Velocity  of  the  Inner  Walls  of  a Liner  Shaped  Charge  Liner," 
Ballistic  Research  Laboratory  Memorandum  Report  No.  1696 , Aberdeen 
Proving  Ground3  MD3  Sep  1965.  (AD  #478326) 

&A.  Kiwan  and  H.  Wisniewski,  "Theory  and  Computations  of  Collapse  and 
Jet  Velocities  of  Metallic  Shaped  Charge  Liners 3"  Ballistic  Research 
Laboratory  Report  No.  1620 3 Aberdeen  Proving  Ground,  MD3  Nov  72. 

(AD  #907161) 
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AXIAL  POSITION  (mm) 


Figure  6.  BRLSC  code  results  of  the  axial  velocity  for  the  liner 
material  along  the  axis  of  symmetry  at  8,  10,  and  13ys 
from  the  43mm  shaped  charge  with  apex  angle  of  45°.  The 
jet  tip  is  on  the  left. 
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Figure  7.  Detailed  analysis  of  calculated  results,  of  the  43mm  shaped 
charge,  along  the  axis  of  symmetry  at'lOys  after  initiation 
of  the  explosive. 
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Interjected  here  are  experimental  observations  regarding  the  continuous 
jet  region,  the  leading  pellet,  for  typical  shaped  charges.  First,  sample 
observations  and  densitometer  measurements  made  from  flash  radiographs 
of  the  jet  in  free  flight  indicate  that  the  density  of  the  continuous 

jet  is  less  than  the  initial  density  of  the  liner  by  about  15  to  20 

percent.7’8  This  density  corresponds  to  a temperature  of  1500°C  which 
is  above  the  melting  point  of  copper.  Thus,  the  liquid  state 
of  the  continuous  jet  is  not  ruled  out,  but  these  measurements 

do  not  take  into  account  the  possibility  of  voids  along  the 

jet's  axis  of  symmetry  which  would  result  in  lower  density  readings. 

Jamet  claims  that  the  jet  breaks  up  in  the  solid  state. 

Second,  radiometry  measurements  of  the  jet's  temperature  were  observed  to 
be  less  than  the  melting  point  of  the  liner  material.9  These  measure- 
ments were  taken  from  both  the  side  on  and  the  front  end  view  of  the 
jet  tip  region.  These  experiments,  conducted  at  the  BRL,  indicate 
that  the  continuous  jet  (tip)  is  cooler  than  the  rest.  Other  experi- 
mental evidence10  indicates  that  the  tip  region  of  the  copper  jet 
breaks  up  sooner  than  the  following  parts.  The  breakup  time  is  faster 
for  higher  strength  material11  and  the  cooler  jet  tip  would  be  expected 
to  have  more  strength.  Calculations  show  that  the  variation  in  tem- 
perature, from  the  jet  elongation  process,  is  negligible  between 
jet  formation  and  particulation.1 1 The  preponderant  evidence  indicates 
that  the  continuous  jet  tip  is  solid. 

Figure  8 a,  b,  and  c are  the  calculated  results  from  the  43mm 
shaped  charge  depicting  the  ratio  of  density,  p/pq;  ratio  of  specific 

internal  energy,  I/Iq ; and  jet  velocity  respectively  for  the  jet  tip 

as  functions  of  time  after  the  initiation  of  the  explosive.  Examining 
these  three  relationships  we  observe  the  following:  first,  the  material 
density  ratio  goes  thru  a maximum  compression  with  rap.id  acceleration 
and  increase  in  specific  internal  energy  for  approximately  2ys;  second. 


7 Private  communication  from  R.  Jameson 3 Ballistic  Research  Laboratory 3 
Aberdeen  Proving  Ground 3 MD. 

SF.  Jamet 3 "Mesure  de  la  densit'e  d'un  get  de  charge  creuse  en 
cuivre  par  radiographic  - e ' clair , " Saint-Louis 3 Rapport  - 
Bericht  R 101/763  9.1.1976. 

3W.  G.  Von  Holle  and  J.  J.  Trimble 3 "Temperature  Measurement  of  Copper 
and  Eutectic  Metal  Shaped  Charge  Jets 3 " Ballistic  Research  Laboratory 
Report  No.  2004 3 Aberdeen  Proving  Ground 3 MD3  August  1977.  (AD  #B021338L) 

l010P.  C.  Chou  and  J%.  Carleone3  "Calculation  of  Shaped  Charge  Jet 
Strain3  Radius  and.  Breakup  Time3"  BRL  Contract  Report  CR2463  July 
197 5,  prepared  by  Dyna  East  Corp.3  Wynnewood3  Pennsylvania. 

(AD#B0072406) . 

1 1i?.  R.  Karpp  and  J.  Simon , "An  Estimate  of  the  Strength  of  a Copper 
Shaped  Charge  Jet  and  the  Effect  of  Strength  on  the  Breakup  of  a 
Stretching  Jet3  " BRL  Report  No.  189S3  June  1976.  (AD#B012141L). 
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1.5 


Figure  8.  Detailed  analysis  of  the  43mm  shaped  charge  at  the  tip  of  the  jet 
giving  three  functional  relationships  of  time  after  the  initiation 
of  the  explosive:  a)  density  ratio*  P /p Q ; b)  specific  internal 

energy  ratio  I/Iq;  c]  axial  velocity,  Vj . 
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the  material  density  ratio  levels  off  to  a constant  with  a gradual 
increase  in  specific  internal  energy  and  velocity  while  equilibrating 
the  inverse  velocity  gradient  for  a period  of  approximately  5ps;  and 
finally,  steady  flow. 


V.  105MM  SHAPED  CHARGE  ANALYSIS 

The  105mm,  unconfined  shaped  charge  data12  is  a valuable  source 
for  comparing  calculated  results.  Figure  9 a and  b are  the  line  draw- 
ings of  the  test  charge  and  the  initial  configuration  used  in  the 
BRLSC  solution  respectively.  Several  reports  have  been  written 
regarding  these  data  and  various  calculational  techniques13,14.  This 
report  is  the  second  effort  to  compare  an  Eulerian,  two-dimensional, 
finite-difference,  hydrodynamic  calculational  technique  against  these 
data.  Comparisons  between  all  currently  existing  calculated  data 
and  experimental  results  will  be  shown,  as  well  as,  calculated 
equations-of-state  variations  in  a parametric  study  of  this  shaped 
charge. 

Figure  10  illustrates  the  classical  collapse  and  jet  velocities  as 
functions  of  relative  distance  from  the  apex  of  the  cone.  The 
respective  final  jet  tip  velocities  of  experiment  and  calculation  are 
7.0mm/ps  and  6.8mm/]js.  This  demonstrates  the  major  focal  point  of 
comparison  between  experimental  and  calculated  results. 

The  procedure  used  in  obtaining  the  calculated  collapse  and  jet 
velocity  data  will  be  presented.  First,  the  collapse  velocity  as  a 
function  of  initial  position  along  the  liner,  was  obtained  by  monitor- 
ing the  center  line  tracer  at  each  time  increment.  As  is  illustrated 
in  Figure  11,  the  center  line  tracer  is  the  third  tracer,  from  five 
tracers,  denoted,  as  T^  thru  TO^,  initially  spaced  in  even  increments 

across  the  liner.  The  value  of  the  collapse  velocity  depicted  in 
Figure  10  is  taken  when  this  monitored  tracer  reaches  its  maximum  as 
a function  of  time.  Then  the  jet  velocity,  as  a function  of  initial 
position  along  the  liner,  was  obtained  by  monitoring  the  inside 
tracers,  (for  example,  T^  and  T^  in  Figure  11),  at  each  time 


12P.  E.  Allison  and  R.  Vitaliy  nAn  application  of  the  Jet-Formation 
Theory  to  a 105mm  Shaped  Charge,  " Ballistic  Research  Laboratories 
Report  No.  1165,  Aberdeen  Proving  Ground,  MI),  March  1961 
(ADU77  458). 

13M.  L.  Gittings,  R.  T.  Sedgwick,  and  J.  M.  Walsh,  " Numerical  Analysis 
of  Jet  formation  from  Lined  Shaped  Charges,  " Ballistics  Research 
Laboratories  Contract  Report  No.  51,  Aberdeen  Proving  Ground,  MD, 
August  1971.  (AD  #889631L) 

14J.  T . Harrison  and  R.  R.  Karpp,  nTerminal  Ballistic  Application  of 
Hydrodynamic  Computer  Code  Calculations,  " Ballistic  Research 
Laboratory  Report  No.  1984,  Aberdeen  Proving  Ground,  MD,  April  1977. 
(AD  #A041065) 
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Composition  B 


Copper  Liner 


DIMENSIONS 
A = 83  mm 


Figure  9a.  The  105mm  unconfined 
test  charge  used  at  the  Ballistic 
Research  Laboratory.  Pertinent 
dimensions  are  shown  in  the  table. 
This  figure  is  from  Ref. 12. 


Figure  9b.  Initial  configuration 
used  in  BRLSC  solution  of  the 
BRL-105mm  unconfined  test  charge. 
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increment.  Therefore,  the  calculated  data  contained  on  Figure  10  is 
the  inside  of  the  liner  tracer's  velocity  distribution  at  the  final 
computed  time  increment,  which  was  30ys  after  the  initiation  of  the 
explosive.  The  calculation  was  terminated  at  that  time  increment 
due  to  initial  grid  boundary  constraints.  At  this  point,  however, 
only  approximately  45%  of  the  initial  liner  had  collapsed  onto  the 
axis  of  symmetry.  The  data  is  presented  as  it  now  stands.  Collapse 
velocities  were  obtained  from  the  entire  liner  since  the  detonation 
wave  has  traversed  the  liner  at  this  time.  Enough  data  are  available 
to  indicate  that  the  calculated  lead  pellet  or  jet  tip  is  projected 
from  the  first  18%  of  cone's  distance  from  the  apex,  where  as, 
experimental  results  show  that  the  lead  pellet  emerged  from  the  first 
39%  of  the  cone's  distance  from  the  apex.  This  difference  is  45% 
of  the  experimental  result.  This  is  illustrated  on  Figure  10  as  a 
constant  jet  velocity  at  the  apex  of  the  liner. 

As  another  measure,  we  are  able  to  compare  the  experimental  and 
computed  radial  component  of  collapse  velocity  as  a function  of  time 
after  initiation  of  the  explosive.  This  is  illustrated  on  Figure  12 
as  the  experimental,  dashed  lines,  and  computed,  solid  line.  The 
experimental  data  were  obtained  by  Randers-Pehrson15 . These  data  points 
are  on  the  inside  surface  of  the  liner  at  various  initial  radii 
positions  monitored  as  a function  of  time  after  the  initiation  of  the 
explosive.  The  computer  code  calculations  are  obtained  by  monitoring 
the  inside  passive  tracer,  for  example  T^  and  T^  on  Figure  11,  as  a 
function  of  time  after  the  initiation  of  the  explosive  from  which  the 
radial  component  of  the  collapse  velocity  was  extracted. 

From  the  two  Figures  10  and  12  we  conclude  that  good  agreement 
exists  between  the  experimental  and  calculated  collapse  velocity 
distributions  and  that  the  jet  velocity  distributions  are  of  the  same 
general  shape.  IVe  further  conclude  that  good  agreement  was  obtained  between 
the  experimental  and  calculated  jet  tip  velocity  with  only  a 3%  varia- 
tion but  the  calculated  lead  pellet's  position  from  the  apex  of  the 
cone  is  45%  of  the  experimental  result. 


VI.  PARAMETRIC  STUDY  OF  THE  105MM  SHAPED  CHARGE  VARYING  THE 

EQUATIONS-OF-STATE 

The  last  study  contained  in  this  report  is  a parametric  analysis 
of  the  105mm,  unconfined  shaped  charge  with  respect  to  the  equation-of- 
state  used  in  this  calculational  technique.  Equations-of-state  for 
both  the  metal  liner  and  the  explosive  filler  were  varied  parametrically. 


1 ^Private  Communi.oati.on  from  G.  Randers  -Pe  hr son 3 Pioatinny  Arsenal 3 
New  Jersey. 
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(S7//WUU)  A1I0013A 


Figure  10.  Comparison  between  experimental  data  from  reference  12  and 
calculated  collapse  and  jet  velocity  distributions  from 
the  105mm  shaped  charge. 
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Figure  11.  Section  of  copper  liner  showing  the  string  of  passive  tracer 
particles  employed  for  predicting  collapse  and  jel  velocity 
results . 


RADIAL  VELOCITY  (mni//is) 


2,0 


Figure  12.  Comparison  between  experimental  data  from  reference  15 
and  calculations  of  the  radial  collapse  velocity  from 
the  105mm  shaped  charge. 
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For  clarity,  the  forms  and  constants  of  these  state  equations  are  pre- 
sented as  appendices  to  this  report.  The  criteria  utilized  for  the 
comparison  are  correlation  of  the  jet  tip  velocity  to  the  observed  value 
of  7.0mm/ys  and  the  flow  field’s  tracer  particle  outlines  at  20ps  after 
the  initiation  of  the  explosive. 

Tabulated  in  Table  1 are  both  the  variable  names  associated  with  each 
equation-of-state  are  their  respective  jet  tip  velocity  results. 

Table  1.  Equations  of  State  and  Jet  Tip  Velocities 


Variation 

Equation 

of  State  Name 

Jet  Tip 

Number 

H.E. 

Metal  Liner 

Velocity  (mm/ps 

1 

Y LA W 

HEMP 

6.2 

2 

Y LAW 

Tillotson 

7.3 

3 

JWL 

HEMP 

6,0 

4 

JWL 

Tillotson 

6.8 

Figures  13-16  are  the  flow  field  tracer  particle  outlines  of  boundaries 
corresponding  to  numbers  1-4  respectively  in  Table  1.  Based  upon  the 
criteria  for  the  comparison,  variation  number  4 was  chosen  as  the  set 
of  state  equations  to  represent  the  best  approximation  of  the  physical 
phenomenon.  This  set  of  state  equations  was  utilized  in  the  previous 
study,  i.e.  the  43mm  shaped  charge  study.  All  plots  resulting 
from  the  calculation  of  variation  number  4 will  be  presented.  The  first 
series  of  plots.  Figures  17-20,  are  the  velocity  fields  at  times  of 
0,  5,  10,  and  20ps.  The  final  series  of  plots.  Figures  21-24  are  the 
two-dimensional  pressure  fields  at  the  same  respective  times  as  above. 

Having  shown  that  variation  number  4 gave  the  best  prediction  of  jet 
tip  velocity,  we  will  compare  it  to  Variation  number  3,  the  worst  tip 
velocity  prediction,  in  more  detail.  The  data  for  variation  number  3 are 
taken  from  Figure  25  which  shows  the  results  of  assuming  a HEMP  equation  of 
state  where  p=f(p).  For  variation  number  4,  Figure  26  shows  the  results 
for  a Tillotson  equation  of  state  where  p=f(p,  I). 

Table  2 summarizes  a comparison  of  the  jet  tip  properties  for 
variations  number  3 and  4,  giving  the  maximum  and  minimum  values.  It  will 
illustrate  the  differences  in  the  computed  results  due  to  the  state 
equations. 
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Table  2. 

Detailed  Analysis 

at  the  Jet  Tip 
Values 

Giving  Minimum 

and  Maximum 

Variation 

P/P 

o 

Velocity  (mm/gs) 

P/PQ 

I/:o 

Number 

(a  minimum) 

(a  maximum) 

(a  minimum) 

(a  maximum) 

3 

P=f(p)=0 

6.0 

1.0 

1.45 

4 

P=f(p,I)=0 

6.8 

.48 

1.45 

This  table  illustrates  that  when  specific  internal  energy  is  used  in 
the  solution  of  the  pressure  equation,  P=f(p,I),  the  jet  tip  is  48%  of 
the  initial  liner  density.  However,  the  jet  tip  is  traveling  at  a 
velocity  which  is  97%  of  the  observed  value.  This  table  also  illustrates 
that  by  eliminating  the  calculated  specific  internal  energy,  i.e.  the 
Hugoniot  relationship,  will  result  in  a density  that  is  equal  to  the 
initial  liner  density,  but,  the  jet  tip  velocity  is  86%  of  the  observed 
value.  As  an  aside,  we  note  that  the  value  of  the  specific  internal 
energy  ratio  at  the  jet  tip  for  these  two  calculations  as  well  as  for 
the  other  calculations  in  the  parameteric  study  are  the  same.  This  value 
is  10.5  times  greater  than  I /IQ. 


VII.  CONCLUSION  AND  SUMMARY 

These  comparisons  show  an  excellent  agreement  between  calculated 
and  experimental  results  during  liner  collapse.  Also,  calculated  jet 
and  liner  collapse  velocities  compare  well  with  experiments.  But,  when 
different  forms  of  equations-of-state  were  applied  in  the  BRLSC  code, 
variations  occur  in  both  the  jet  density  and  velocity.  The  specific 
internal  energy  seems  to  be  the  source  of  the  variations  in  the  calcula- 
tions. If  these  hydrodynamic  computer  codes  are  ever  going  to  be  able 
to  model  the  entire  physical  process  with  detail  and  accuracy,  then 
additional  groundwork  must  be  accomplished.  We  suggest  investigating 
the  following: 

a.  Check  the  equations-of-state, 

b.  determine  if  the  numerical  scheme's  zone  size  is  responsible 
for  the  exaggerated  over-thermalization  of  specific  internal  energy, 

c.  determine  whether  or  not  a new  numerical  scheme,  i.e.,  a second 
order  or  possibly  a simple  modification  to  the  existing  scheme,  will 
bring  better  results. 

As  mentioned  earlier,  these  determinations  must  be  made  before 
we  can  proceed  with  more  sophisticated  investigations  such  as  strength 
effects  and  jet  breakup. 
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Figure  13.  Computed  flow  field's  tracer  particle  outlines  of  results 
from  case  number  1 on  Table  1.  Results  are  plotted  at 

20ys  after  initiation  of  the  explosive. 
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Figure  14.  Computed  flow  field's  tracer  particle  outlines  of  results 

from  case  number  2 on  Table  1.  Results  are  plotted  at  20ys 
after  initiation  of  the  explosive. 
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Figure  15.  Computed  flow  field's  tracer  particle  outlines  of  results 

from  case  number  3 on  Table  1.  Results  are  plotted  at  20ps 
after  initiation  of  the  explosive. 
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Figure  16.  Computed  flow  field's  tracer  particle  outlines  of  results 

from  case  number  4 on  Table  1.  Results  are  plotted  at  20ys 
after  initiation  of  the  explosive. 
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Figure  17.  Velocity  Field  in  the  Eulerian  grid  at  time  = 0 ys.  Results  from  number  4 on  Table 
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Figure  18.  Velocity  field  on  the  Eulerian  grid  at  time  = 5ps.  Results  from  run  number  4 on  Table 
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Figure  19.  Velocity  field  on  the  Eulerian  grid  at  time  = 10ys.  Results  from  run  number  4 on  Table 
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Figure  20.  Velocity  field  on  the  Eulerian  grid  at  time  = 20ys.  Results  from  run  number  4 on  Table 
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Figure  21.  Three  dimensional  pressure  field  on  the  Eulerian  grid  at  time  = Ops.  Results  from 

run  number  4 on  Table  1. 
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Figure  22.  Three  dimensional  pressure  field  on  the  Eulerian  grid  at  time  = Bps. 

Results  from  run  numhpr  4 nn  Table  1. 
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Figure  24.  Three  dimensional  pressure  field  on  the  Eulerian  grid  at  time  = 20ys.  Results  from 

run  number  4 on  Table  1. 


a. 


c. 


Figure  25.  Detailed  analysis  of  calculated  results  of  the  105mm 

shaped  charge  (see  number  3,  on  Table  I),  along  the  axis 
of  symmetry  at  20ys  after  initiation  of  the  explosive. 
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Figure  26.  Detailed  analysis  of  calculated  results  of  the  105mm 
shaped  charge  (see  number  4,  on  Table  I),  along  the 
axis  of  symmetry  at  20p  after  initiation  of  the  explosive 
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APPENDIX  A 


I . Equations  of  State  for  Copper  Liner 

There  are  two  forms  of  the  equation  of  state  used  for  the  copper 
liner.  The  first,  used  for  determining  pressures  and  constant  energy 
compressibilities  in  the  liner,  is  that  due  to  J.  H.  Tillotson16. 

Only  the  condensed  form  of  the  equation  of  state  is  included  since  the 
liner  material  should  not  be  shocked  enough  to  cause  significant 
vaporization  of  material.  The  equation  of  state  formulations  are, 

for  n > 1 - Jt  : 


P = 


a + 


+ 1 


Ip  + Ay  + By  , 


and  when  n < 1 - ^ 

ZD 


P 


b 

— 

I 

+ 1 

T 2 

*o  n 

If> 


In  these  equations  P,  I,  and  p are  pressure,  specific  internal  energy, 

and  density  respectively,  n = — = y + 1,  and  the  constants 

po 

for  the  copper  liner  are  as  follows: 

a = .5 

b - 1.5 

A = 1.39  Mbar 

B = 1.1  Mbar 

I = .325  x 103  joules/g 

pQ  = 8.9  Mg/m3 


16J.  H.  Tillotson 3 "Metallic  Equations  of  State  for  Hyper-velocity 
Impacts "j  General  Atomic  Report  GS-3216 3 July  1961. 
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The  second  is  taken  from  the  HEMP  code  report17.  The  equation-of-state 
formulations  are, 

for  h > 1 - 4s-  : 


P 


= Ay  + By 


2 


J 


and  when  n < 1 - 


A_  . 
2B  ' 


The  constants  for  the  copper  liner  are  as  follows: 

A = 1.19  Mbar 
B = 4.35  Mbar 
PQ  = 8.9  Mg/m^ 


^Unpublished  HEMP  User's  Code  Report  from  M.  Wilkens3  Lawrence 
Radiation  Laboratory. 
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APPENDIX  B 


II.  Equation  of  State  for  Composition  B Explosive 

There  are  two  forms  of  the  equation  of  state  used  for  the  grade  A, 
composition  B explosive.  The  first  is  a modified  gamma-law  equation 
of  state  of  the  explosive. 


= (Y  - l)p  MAX[I 


(I  - C • I )] 


min 


In  this  equation,  P,  I,  and  p are  the  pressure,  specific  internal 
energy,  and  mass  density  respectively.  The  constants  for  composition  B 
explosive  for  y,  IQ,  and  p and  2.71,  5.034  x 103  joules/g,  and  1.717 
Mg/m3  respectively.  1,^  is  an  input  parameter  (usually  about 
10"1  joules/g).  The  variable  C is  a flag  ranging  between  0 and  1 which 
determines  whether  a cell  is  behind  (C  = 0) , intercepted  by  (0  < C < 1) 
or  ahead  (C  = 1)  of  the  detonation  front.  For  details  see  report 
by  M.  L.  Gittingsl. 

The  second  uses  the  Jones-Wilkins-Lee18  (JWL)  equation  of  state  for 
the  grade  A,  composition  B explosive. 


P = 


(‘  ' “O' 


+ ap  MAX 


ram 


Bp, 


(I  - C 


where:  A = 5.2423  Mbar 

B = .076783  Mbar 

a = .34 

a = 4.2 

B = 1.1 

I = 4.95  x 103  joules/g 

pn  = 1.717  Mg/m3 


The  variable  C and  constant  IQ  are  used  to  determine  the  cells 
intercepted  by  the  detonation  front  as  explained  above. 


18i?.  Lee,  M.  Finger , and  W.  Collins , "JWL  Equation  of  State  Coefficients 
for  High  Explosives 3 " Lawrence  Livermore  Laboratory , UCID-16189 , 
January  196 3. 
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